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ABSTRACT 


A summary of previously proposed receiver strategies for direct-detection laser ranging re- 
ceivers is presented. Computer simulations are used to compare performance of candidate 
implementation strategies in the 1- to 100-photoelectron region. Under the condition of no 
background radiation, the maximum-likelihood (ML) and minimum mean-square error 
(MMSE) estimators were found to give the same performance for both bell-shaped and rec- 
tangular optical-pulse shapes. For signal energies greater than 10 photoelectrons, the root- 
mean-square (rms) range error is shown to decrease as foj- bell-shaped pulses and Q'* 
for rectangular pulses, where Q represents the average pulse energy. Of several receiver 
implementations presented, the matched-filter peak detector was found to be preferable. 

A similar configuration, using a constant-fraction discriminator, exhibited a signal-level 
dependent time bias. Suggestions for future study are also included. 
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A COMPARATIVE STUDY OF OPTIMUM AND SUBOPTIMUM 
DIRECT<DETECTION LASER RANGING RECEIVERS 


James B. Abshire 

Goddard Space Flight Center 
Greenbelt, Maryland 


INTRODUCTION AND SUMMARY 

Since the development of high-power pulsed lasers, there has been a widespread interest in 
the application of these devices to optical radar systems. The development of laser ranging 
systems within the National Aeronautics and Space Administration (NASA), for example, 
began in the early 1960’s and has grown substantially in the intervening years (Reference 1). 
The application of accurate range measurements to Earth-orbiting satellites has produced 
important results in precise orbit and gravity-field determinations and is expected to produce 
in the near future the same unique findings relating to the crustal motions of the Earth. 
Because the potential usefulness of laser ranging data is directly related to its accuracy, there 
is considerable interest in developing optimum laser ranging receivers. 

Several studies have been made of optical radar receivers. An early study (Reference 2) 
compared the performance of direct-detection versus heterodyne receivers and found that, 
for conditions frequently encountered in practice, the direct-detection receiver had superior 
target detection and ranging performance despite the greater complexity of the heterodyne 
receiver. The performance of the heterodyne receiver was found to exceed the simpler 
direct-detection receiver only when background noise was very heavy or when high-velocity 
resolution was required. 

As a result of these findings, several methods have been outlined for implementing the direct- 
detection receiver when the return optical pulse shape is known in advance (References 3 
through 5). In addition, recent work has studied direct-detection receivers for use when the 
return pulse shape is randomly broadened by the radar target (References 6 and 7) 

In this document, the ranging performance of several proposed direct-detection receivers 
is compared. As necessary background, the minimum mean-square error (MMSE) and 
maximum-likelihood (ML) estimators for optical pulse delay are derived. These estimators 
are shown to be equivalent in the special case of the receiver operating with no background 
illumination and with either bell-shaped or rectangular pulse shapes. The error performance 
of the ML estimator is derived under the assumption of high optical signal to background 


ratios and differentiable optical pulse shapes. The results of computer simulations show the 
performance of the ML receiver operating under low signal energies and with increasingly 
rectangular pulse shapes and represent new results in this area. The simulated performance 
of a modified ML receiver is also presented. 

Finally, the ranging performance of optimum and suboptimum receiver configurations is 
compared, and the results are used to outline possible improvements to current laser-ranging 
receiver design. 

MATHEMATICAL REPRESENTATION OF A POISSON PROCESS 

It has been shown that the beam fluctuations, and therefore the photocount, from a single- 
mode laser beam follow Poisson statistics (Reference 8). Furthermore, when noise is absent, 
the photocount rate varies as the signal that modulates the beam and, when noise is present, 
as the signal plus noise. 

The noise in laser-ranging receivers is attributable to both detector dark current and back- 
ground radiation, which are also believed to follow Poisson statistics. In all NASA laser ranging 
systems that have been planned or built to date, high-gain photomultipliers have been used 
as detectors. In systems using these detectors, effects of thermal noise on receiver perfor- 
mance are negligible. 

For the theoretical model used in this study, the laser transmitter is treated as a photon 
source whose average output power is modulated by a deterministic waveform (t) of 
limited time duration. In general, the received data will consist of a series of randomly 
distributed points in time, t^, j = 1 ,2, . . ., M, within the time interval, (-T, T). Here, both 
the t.’s and the total number of observed photocounts, M, are random variables. The t2s 
are a realization of a Poisson point process with time-varying expectation X (t). For the 
range estimation problem, X (t) = X^, (t - t) + X^^ , where X^, (t) is known and X^ represents a 
uniform background noise count rate (attributable to background illumination on the photo- 
detector or internally generated detector noise counts). The time interval, (-T, T), represents 
the limits on the prior range uncertainty and is assumed to be equal to the observation interval. 
For the MMSE and ML estimators, it is assumed X^. (t - r) = X^ (t - t) (that is, the average 
received optical pulse shape is equal to that of the transmitter). The receiver is required to 
estimate the time-of-arrival of the optical pulse, t, given observation of the photocounts, t., 
j = 1 ,2, . . ., M within (-T, T). 

This simple model of the optical detection process has at least two limitations. First, if the 
optical pulse was produced by several interacting laser modes or by an incoherent optical 
source, the photoelectron distribution is not necessarily Poisson (Reference 8). Also, the 
model assumes that observation of the individual photoelectric emission times at the detector 
is accurate. In practice, photomultipliers used as detectors in laser-ranging systems have time 
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jitters ranging from 30 to 700 picoseconds (ps) at the single photoelectron level. Such photo- 
multiphers also exhibit random gains on a pulse-by-pulse basis, which can cause uncertainties 
in the number of photocounts that occurred within the resolution time of the detector 
(Reference 9). Because the effects of stochastic detector transit times and gains have not 
been incorporated in more elaborate detector models to date, the photodetector is assumed 
to be ideal in the following analysis. 

To derive the probability density functions of the foregoing Poisson process, it is assumed 
that the following axioms are valid: 

• The probability of one photocount in an infinitesimal interval of width, At, is 
given by P [ 1, At] , where 


Lim P [1, At] = X(t)At (1) 

At-^O 

• The probability of more than one photocount in the time interval, At, approaches 
zero as At ^ 0 

• The number of photocounts in any one interval is independent of those for all 
other disjoint intervals. 

Given these axioms, the probability that the number of photocounts in interval (a, b) equals 
integer K is given by 


P [K, (a, b)] 



where X (t) > 0 is the time-varying expectation of the Poisson process. 
For K = 0 and K = 1, equation 2 reduces to 


P[0,(a.b)] 



( 2 ) 


(3) 
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X (t) dt 


(4) 


b 

P[l,(a,b)] = P[0,(a,b)] . y* 

a 


Consider a realization of the process as shown in figure 1. 


^ At, ■«- ->• Atj ■«- -»■ At^^ <- 

— j 1 — •— j 1 — •—! — • — ••• — • 1 — • — |- 

-T *1 h *M-i Va 

Figure 1. A possible realization of the random-point process. 

Let I tj^ denote the set of M numbers, - T < t , < t^ < ' * ’ < tw ^ 

By first allowing infinitesimal intervals of width, At^, about t^ j = 1,2,.. ., M, the M 
dimensional-joint probability density function, P (|tj^[ , M), can be calculated as follows. 
By axiom 3 and equation 4, the probability that only one point will fall within each of 
these intervals with no counts outside of them is given by 


-j ^ TIME 

T 


P(|t^,},M)-[At, -Atj *... -Atj^] =[P[0,(-T,t, - 1/2 At,)] .X(t,)At,] 

. [P[0,(t, + l/2At,,t2 - l/2At2)j -X(t2)At2] 

• [P[0,(V, + l/2At,„.,,t,„ - l/2AtJ].X(tM)AtM] 

• P[0,(t^, + l/2At^,,T)] 


By equation 3, in the limit as At^ 0 for all j and M > 1, 


P({t,^},M) = exp 



^ M 

X(t)dt nx(t.) 

) j=i 


(5) 
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If no occurrences are observed, M = 0, | tj^ | is empty, and 


P({t„},0) = P[0,(-T,T)] = e-Q 


( 6 ) 


where 


Q = 



(7) 


Figure 2 illustrates the sample space of events for this random-point process. In figure 2, 



Figure 2. Sample space for the random-point process. 

denotes the set of no occurrences in (-T, T); , the set of one occurrence in (-T, T); 

, the set of all possible sets of M events {t^ } ; and the sure event (i.e., some number 
of points occurred at any instants in the interval (-T, T)). 

As a special case, consider the process to be stationary with constant mean X (t) = K. Then 
all sample sequences of M occurrences within the interval, (-T, T), have the same likelihood, 
with joint probability density 


P({tM}.M) = (K)Mexp|- 2KT} 


( 8 ) 
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DERIVATION OF OPTIMUM ESTIMATORS 


The mathematical forms of the MMSE and ML estimators are derived in this section, follow- 
ing methods previously outlined (References 3 and 4). These forms are shown to be equiv- 
alent for the case of Gaussian optical-pulse shapes and no background noise. 


Minimum Mean-Square Error Estimators 

For a given set of observed data, it is well-known that the conditional mean estimator 
minimizes the mean-square error of the estimate. However, to find the conditional mean 
estimator, the value of the optical-pulse delay, t, must also be assumed to be a random 
variable. In addition, the prior density function of the delay, (r), must be known over 
the observation interval, (-T, T). This requirement can be severe for ranging receivers because 
often little is known about the statistical fluctuations in the range of the target. For example, 
in satellite laser ranging, the major purpose of the experiment is to gain just this statistical 
knowledge concerning the ranging station to satellite distance (Reference 1). Nonetheless, 
it is still instructive to derive the MMSE estimator, because its functional form can be given 
explicitly for cases that cannot be solved by the ML estimator, as the next section will show. 
For simplicity in the following derivations, the prior range density, P (t), is assumed to be 

T 

uniform over (-T, T), and the MMSE estimator is solved for rectangular, exponential, and 
Gaussian optical-pulse shapes. 

As before, let r denote the unknown range delay, and let P^ (t) denote its prior probability 
density. Also, let t, j - 1,2, . . ., M denote the ordered instants of photoelectron emissions, 
and let P (|t^j [ It) denote the conditional density of | , given t. The MMSE estimator 
f is the conditional mean of the range, given the observed data, and can be written as 


where 


- f rP(|t„(lr)P,Md, 


P(|tM},M) = / 


(9) 


( 10 ) 


and R = the observation interval. 

T 
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Assuming Poisson statistics and X (t) for the optical pulse shape and using equations 5 
through 7, 


P(K(lr). 


M 

e-^ n X(t. ~ T) ,u>\ 

j=l 


,M = 0 


Substituting equations 10 and 1 1 into equation 9 yields 


( 11 ) 


r = 


J' rj^n X(tj-T)Jp^(T)dT 

/ " M 

n X(t.-r) 

R. Lj=l J 


(r) dr 


(12) 


Note that, by this definition, f is undefined when no photoelectron emissions are observed. 
The equation 12 will be solved for rectangular, exponential, and Gaussian pulse shapes in 
the background free case. 


Rectangular Pulse Shapes 

Let the pulse shape be rectangular of width 2D, centered about zero. Then 


and 


X(t) 


i Q/2D = Xs 

( ° 


, It KD 


, otherwise 


(13) 


•■({Ui') 



, It.^Tl <D,fora]l j<M 
,lt.-T|>D 


(14) 
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If the prior density of r is uniform, 


j(2T)-i JtKT 

P,(r) = 

V 0 , otherwise 


(15) 


Therefore, 


■’(tU.M) 


y e-° Xs“ dr - Xs“ 
R 2T 2T 

T 



(16) 


Figure 3 illustrates the method used to find R^, the range of the unknown delay value given 
1 1^ I . Figure 3 shows that r has the range, [tj^ " + D] • 


/ 


OPTICAL PULSE 
ENVELOPE 



1 


1 1 
1 

T 

max 

1 



OBSERVED . 
DATA 


• ♦♦♦ 




TIME 


mn 


Figure 3. Prior range uncertainty limits given M photocount observations. 
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Using these limits for R^, 




e-Q \ 
2T 


M 


(2D + - t„) 


- tj<2D 


(17) 


, otherwise 


ForM> 1 andt^, -tj<2D, 


/ ''■(I'mIi 


='«X“ 

r)P^(r)dr = 




(2D + tj - tjj) 


(18) 


Substituting equations 1 7 and 1 8 into equation 9 yields, for rectangular pulses, 




rect 


(19) 


which is valid for M > 1. 

Exponential Pulse Shapes 

Let the pulse shape be a decaying exponential with its origin at t = 0. Then, 


X(t) = 


^ exp (- t/D) , t > 0 


, otherwise 


Assuming that t > 0 and M > 1, 


( 20 ) 


M M 

P(^tM}lr) = e-Q n X(t. - T) = e-9(Q/D)“ Yl e'V^'^ 
j=i j=i 


= e-^ (Q/D)“ exp 




( 21 ) 
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where 


M 

t = (M)-» 2^ t. 

j=l 


Substituting equations 20 and 15 into equation 10 yields 






Mr/D 


dr 


( 22 ) 


From analysis similar to that shown in figure 3, the range ofr is found to be [-T, tj ]. 
Therefore, equation 22 can be evaluated as 





-Mt/D 

e 




(23) 


where 


a = (Q/D)“ exp(- Q - M t/D) 


Similarly, the numerator of equation 12 can be evaluated as 


ot 



^eMr/D 


dr = a 


Mt/D 

e 



\ -MT/D / DT \ 

^ ^ \ M ’ / 


(24) 


Therefore, given M > 1, the optimum estimator for exponential pulse shapes is 


Mt,/D -MT/D 

(tj-D/M)e ^ +(T+D/M)e 


Mt/D -MT/D 

e ^ - e 


(25) 
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Note that, as for most cases of interest, D « T and equation 25 simplifies to 


■ (*1 - D/M) 


( 26 ) 


Bell -Shaped (Gaussian) Pulse Shapes 

Let the optical pulse shape be Gaussian, centered at t = 0 with variance = . Then, 



(27) 


Assuming that M > 1, 



Substituting equations 28 and 15 into equation 10 yields 




f-i— 


M 


M 


y|-) 


dr 


j=i 


(29) 


= ^ 




dr 
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where 


P = — (-p^l ,andR^ = [-T,T] 
2T\v'27rD/ 


Expanding the quadratic term in equation 29, 





T 


2 T c + M r^) 


d T 


(30) 


where 


M M 

b = ^ and c = ^ t . 

j=i j=i 

Because D « T for the case of interest, the area of the Gaussian pulse outside range (-T, T) 
is negligible. Therefore, 





(b - 2tc + Mt^) 


dr 


which can be evaluated to yield (Reference 1 0) 


P(K}.M) = ^ 



(31) 


( 32 ) 


In a similar fashion, the numerator of equation 12 becomes 



which can be evaluated as (Reference 11) 


cDexp 


- bM 
2D^ M 


(33) 


For M> I, the ratio of equation 33 to equation 32 yields the optimum estimator for a 
Gaussian pulse. 


^GAUSS 



(34) 


This estimator is simply the average of the observed emission times. 


Maximum-Likelihood Estimator 

When the prior distribution of ranges is completely unknown, an optimum estimator can 
still be formed by directly optimizing P (|t^{^ , M) for the observed photo-occurrence times, 
{t|^ f . In this section, the ML estimator is derived, and its form is compared to the MMSE 
estimators derived in the previous section. 

In general, let the time-varying expectation of the Poisson process depend on the parameter 
a so that 


X(t,a) = S(t,a) + X„ 


(35) 


where X„ represents a constant background count rate, and S (t, a) is the optical pulse shape. 


which is assumed to be known in both t and a. 
Using this form for X (t), equation 7 becomes 

T 

Q(«)= f 

-T 


Here, a represents a general parameter. 


X(t,a)dt (36) 
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and represents the average number of photocounts observed in interval (-T, T). The likeli- 
hood function of the M observed photocounts can be written as 

M 

P({tM}.M) = e-Q(“) n [S(t.,a) + XJ (37) 


The ML estimate of a is that value of a denoted as a, which maximizes equation 37. Because 
the logarithm is a monotone nondecreasing function of its parameter, 


(38) 


In general, a search procedure is required for finding a. However, if S (t, a) is piecewise 
differentiable, a must satisfy 


a = Max^ I S (V 


_a 

da 


M 


j=l 


£n [S (t., a) + XJ 


Q(a) 


= 0 


which implies that 


_d_ 

da 


Q(a) 


E 

izz 1 



S(t,a))/[S (t., a) 



(39) 


The ML estimate of a must satisfy equation 39 for the general estimation problem. 

For ranging applications, the parameter, a, is delay time, r, and the optical pulse shape can 
be written as X^ (t - t) as defined earlier. For this case, equation 38 becomes 


! m 

^Cn[X^(t.- r) + X^] - Q(t) 


(40) 
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For most cases of interest, the signal pulse is well contained within the prior range uncertainty, 
which implies that the average photocount rate, Q, does not vary with delay time, t. There- 
fore, if X (t) is once differentiable, the ML estimate of range delay, r, must satisfy 




(41) 


assuming the estimator has only one maximum. 

For more general pulse shapes, the ML estimate, f, can be found by maximizing the log 
likelihood ratio, 


M 

F(r) = «n[X^(t.-r) + Xj (42) 

j=i 

with respect to t. This distinction between differentiable and nondifferentiable pulse 
envelopes is important, because the error expression for the ML estimator can be found 
analytically for differentiable pulse envelopes only. For nondifferentiable pulses, other 
techniques are necessary for estimating the ranging receiver performance. 


Comparison of Optimum Estimators 

Table 1 summarizes the form of the optimum estimators of range delay for estimating 
MMSE and ML derived previously. In this section, the forms of these expressions will be 
examined with regard to their hardware implementation. 

To directly implement the MMSE estimator expression for the Gaussian pulse, each photo- 
electron emission time must be observed and recorded, and the arithmetic mean must be 
computed. For nanosecond pulse-width lasers, this processing of the individual photo- 
occurrence times is only possible with streak cameras or swept-image converters. Although 
such high-time resolution devices are now being used in laboratory applications, their relative 
complexity places this form of estimation at a disadvantage relative to other forms. In con- 
trast, the expression for ML estimators for differentiable pulses can be readily implemented 
in two forms. By constructing the time functions. 


q(t) = X,(t)/[\(t) + XJ 


(43) 
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Table 1 

Optimum Estimator Forms for Differentiable 
and Nondifferentiable Pulse Envelopes 


Pulse Type 

MMSE 


Differentiable 

Gaussian: 

General: 

pulse envelopes 

M 

j=i 

M 

j=i 



f = Max,. 

Nondiff eren liable 

Rectangular: 

General: 

pulse envelopes 

^ )/2 

f = Max 

T 


Exponential: 



f = tj - D/M 



ML 


= 0 


or: 


M 

= Max^ y^|en[X3(tj-T) + X„] }■ 

j=l 


M 


j=l 


(«n[X3(tj-r) + Xj| 


and 


V (t) = Cn [X^ (t) + X^] 

either a correlation or a matched filter receiver can be constructed as shown in figures 4 and 5. 

The instant when the output passes through zero (for q (t)) or a maximum (for v (t)) yields 
the estimate, f. The start pulse to the time interval unit is usually derived from the trans- 
mitted laser pulse under extremely high signal conditions. Therefore, the range error intro- 
duced by the start channel is negligible in current system designs and will not be considered 
further here. 

When such a system is operating with low signal-to-background ratios in the return signal 
path, several maxima or zero crossings can occur, yielding ambiguities in the location of the 
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q<t) 

(or v(t)) 





Figure 4. Correlation type ML receiver. 


PHOTONS 


MATCHED 

FILTER 




h(t) = q(L-t) 


1 

or 

h(t) = v(L-t) 


PHOTON 

DETECTOR 


ZERO 
(OR PEAK) 
DETECTOR 


TIU 

STOP START - 


L = DURATION OF SIGNAL ENVELOPE 


Figure 5. Matched-filter ML receiver. 

range estimate. Under these conditions, the maxima that appears in the neighborhood of 
the densest photo-occurrences should be selected. An error in selecting this neighborhood 
can cause a ranging error much larger than an error in locating the occurrence time of the 
zero crossing or maximum. An expression for the probability of this event has been previously 
derived (Reference 3). Because of the relatively high signal-to-background count-rate ratio 
in current NASA ranging systems, this error should occur infrequently. However, for opera- 
ting under less favorable ranging conditions, the increased probability of this event must be 
accounted for in the receiver design. 

For rectangular and exponential pulse envelopes, the expressions for the optimum estimators 
are also directly implementable as shown in figures 6 and 7. 



INVERTER 


Figure 6. MMSE receiver for rectangular pulse envelopes. 
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For rectangular pulse envelopes, a low-pass filter is used following the photodetector. The 
impulse response, h (t), of this filter should be long enough to prevent its output voltage 
from returning to zero between successive photoevents within the duration of the optical 
pulse. The first discriminator will trigger on the leading edge of the filter output, which 
occurs at time tj . The trailing edge of the filter output will trigger the second discriminator 
at time tj^. The threshold levels of both discriminators should be set high enough so that 
they do not trigger on the widely spaced single-photoelectron background counts. The out- 
put at time t^ will stop the range-time interval counter, whereas the time between the first 
and last photoemissions will be measxired and digitized by a time-to-pulse-height converter 
followed by the analog-to-digital converter. With these two measured quantities, a small 
minicomputer or microprocessor can easily compute the correction to tj and form the 
optimum estimator, t. 


/Od. 


ADC 





LOW- 

PASS 

FILTER 


LEADING 

EDGE 

DISC. 


PHOTON 

DETECTOR 


START 


? = t, - D/M 


COMPUTER 


Figure 7. MMSE receiver for exponential pulse envelopes. 


For exponential pulses, only one discriminator is used following the low-pass filter. Although 
the operation of this configuration is similar to that mentioned previously, the total number 
of photocounts, rather than t^^ , is required as the second measurement. By passing the low- 
pass filter output through an integrator, its output at the end of the observation interval is 
directly proportional to M. This output is then digitized, and, as before, the computer can 
form the optimum estimator, f. 

A comparison of the four possible receiver implementations shows that, for receiver sim- 
plicity, one of the ML receivers should be chosen. Furthermore, it is currently more practical 
to implement the matched-filter type receiver than the correlation type, which requires a 
high-speed randomly triggered sampler unit. As was mentioned earlier, the matched filter 
can be constructed for one of two possible waveforms. For maximum flexibility with regard 
to input optical-pulse envelopes, it is evident from table 1 that the filter should be matched 
to V (t) and followed by a peak detector. This configuration is shown in figure 8, and its 
ranging performance is discussed in the following section. 
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A variation of this receiver, using a 50-percent risetime or constant-fraction discriminator 
in place of the peak detector is discussed in the section on “Performance of Suboptimum 
Receiver Strategies.” 


PHOTONS K 


MATCHED 

FILTER 




m 


IMPULSE 
RESPONSE: 
h(t) = v(L-t) 


PEAK 

f 

START 

Til 1 

I 



DETECTOR 

STOP 

1 lU 


PHOTON 

DETECTOR L = DURATION OF SIGNAL ENVELOPE 


Figure 8. Implementable optimum ranging receiver configuration. 


OPTIMUM AND SUBOPTIMUM ESTIMATOR PERFORMANCE 

Because hardware considerations indicate that the ML form is the preferred optimum 
receiver, attention is first focused on deriving the error performance of this estimator. 

Although the following treatment closely resembles that of Bar-David (Reference 3) and 
gives only asymptotically valid results for differentiable pulse envelopes, it is nonetheless 
useful for a wide variety of pulse shapes encountered in practice. However, because the 
chosen hardware implementation is valid for general pulse shapes, results are cited that show 
the performance of this receiver used with rectangular and exponential pulses. Next, the 
results of a computer simulation of receiver performance are described. These results show 
how receiver performance changes when going from bell-shaped to rectangular optical pulse 
shapes, and when using leading-edge rather than peak detection following the matched filter. 
Finally, the ranging performance of these configurations is compared to that of a simple 
leading-edge receiver, which might be used when the optical-pulse shape is randomly broadened 
by the ranging target 


Variance of the ML Estimator 

The variance of the estimate of delay, f , can be approximated by expanding the log-likelihood 
ratio, F (r), given in equation 42 at the true value of delay, . If the first two derivatives of 
F (r) exist, equation 42 can be expanded in a Taylor series as 

• • 

F (T) = F (r„) + F (r„) (r > r„) + ^ (T„)(r - + - (44) 
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Because the ML estimate is the value of t that maximizes F (r), 


9t 


F(r) 


= F(f) = 0 

r=f 


( 45 ) 


Differentiating equation 44 with respect to t, and using equation 45 yields 


F (;> = 0 = F (r„) + (r - r„) F ( t ^) + - 


(46) 


from which it follows that the error in the estimate is 


= ? - To = - F (to)/F (to) 


(47) 


provided F 0. This constraint excludes strictly rectangular pulses from this analysis. 

Thus, the mean value and variance of the delay estimator error are 


E [ej = - E [F(ro)/F(To)] 


(48) 


and 


Var[eJ= E [(F (To)/F (to))"] - [E(e,)]2 


(49) 


where E [ ] denotes the'M-fold expectation over the set of random variables, | tj^ [ . By 
assuming (E [F (To)] )^ » Var [F (Tq)] , F(To) can be replaced by its expected value in 
equations 48 and 49. Under this assumption. 


E[eJ = - E[F(To)]/E[F(ro)] (50) 

and 

Var [e^] = | E [(F (To))^] - 
= Var I F(To) |/E^ 


E"' [F(ro)]|/E2 [F(To)] 


(51) 
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The characteristic functions of F (r^) and F (r^) are 


Char [F ( t ^)] 



= E 


= E 



- ^o)/\ - ^o) + 


and 


Char [F (r„)] 


= E|exp|ix ^ + ^n)]} 

/ 9 . l" 

n exp I ix — \ (t. - 7p)/X^ (t. - Tp) + \„) j. 


• (52) 


(53) 


This can be evaluated by using Bar-David’s results regarding expectation of product functions 
(Reference 3): 


M m — 

n f(t.) = e-Q 2^{[G(T)]^IW.) 

j=l m=0 


g-Q + G(T) 


(54) 


where 


G(T) = 


T 


/ 

-T 


f(t)X(t)dt 


(55) 
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Without loss of generality, the true delay, can be set to zero, and, using equations 54 
and 55, 


Char [F (0)] = exp 



(56) 


Using the welhknown property of characteristic functions. 


E [F(0)] = 


-^Char [F (0)] 
ox 


x = 0 


- i ( 


Char [F (0)] 



x=0 




Char [F (0)] 



X (t) dt 


Because equation 56 is equal to 1 when evaluated at x = 0, 


E[F(0)] = 


1 

/ 

-T 


X (t) dt 


(57) 


Assuming that the pulse is well contained in (-T, T) implies that 


X(±T) = X^(±T) = X^(±T) = 0 


which implies that 


E [F (0)] = 0 


(58) 
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Thus, by equation 50, the ML estimator is unbiased. 

Proceeding in a similar manner with the second derivative of F (r), 


Char [F (0)] = exp { - Q + / X (t) exp i ix F (t)| dt ' 


1 

/ 


( 59 ) 


To find the mean value of F (0), 


[F(0)1 = 


= Char [F (0)] 


lx=0 


-T 


T 

■/[»»- 4#1 


y* F(t)x(t)dt = j |(x(t) - idt 


which simplifies to 


E[F(0)] = 


= -/ 

-T 


(Mi)? 

X(t) 


dt 


(60) 


To complete the calculations for Var (e^), 


Var[F(0)] = 


9x" 


Char [F (0)] 


1 x=0 


= [Char F (0)] | J' [F(t)] ^ X(t) exp ] ix F(t) \ dt j 


x=0 


Using equations 56 and 58, the foregoing expression simplifies to 


Var[F(0)] = 



[X(t)]^ 

X(t) 


dt = - E [F (0)] 
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Finally, by using the definition of X(t), 


Var(e^) = 




-1 


dt 


for the variance of the ML estimator timing error. 


ML Estimator Error for Bell-Shaped Optical Pulses 

As an application of these results, let the pulse shape be a raised cosine, 


\ (t) = 



-D D 

<t< 

2 2 


, otherwise 


(61) 


(62) 


where Q represents the average photoelectron level of the received optical pulse. This pulse 
shape closely approximates a Gaussian pulse shape, whose ±3a points lie at ±D/2, The ML 
estimator variance for a ranging system using such a pulse shape is derived from equation 61 : 


Var(e ) = — 
D 




sin^ (2irt/D) 


1 + cos (iTTtlD) + 


Xn 


dt 


-1 
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where is the average number of background counts per pulsewidth, D. This result is 
well known and has been obtained by others (Reference 5). Note that, if the background 
noise is zero, the root-mean-square (rms) range error has the form 


= D/[27ry^] 


(64) 


where D is the full width of the optical pulse, and Q is the average number of photoelectrons 
per pulse. 


ML Estimator Performance for Increasingly Rectangular Pulse Shapes 

As shown earlier, the performance of the ML estimator cannot be found directly for systems 
using nondifferentiable optical pulse shapes^ For practical systems, this does not present a 
fundamental problem because infinite bandwidth optical sources do not exist; therefore, 
all derivatives of real pulse shapes will be bounded. Nonetheless, because real pulse shapes 
are often more closely approximated by nondifferentiable envelopes than by differentiable 
ones, it is worthwhile to examine the ML estimator error performance of nondifferentiable 
pulse shapes. 

It has been found that under conditions of no background with bell-shaped optical pulses, 
both ML and MMSE estimators asymptotically achieve the Cramer-Rao lower bound.* An 
expression for ML estimator performance was also found: 

Q e-Q “ Q>^ 

= B 7 (65) 

Q 1 - e'^ k ' k! 

k=l 


Var [e^l m > 1] 


which is based on at least one photoelectron occurrence. Here, Up represents the second 
moment of the pulse shape. As Q becomes large, given a raised cosine pulse shape, equation 
65 reduces to equation 63, where X^ is taken as zero. 

To further investigate the ML estimator performance, a computer program was written that 
simulates the behavior of an ideal (i.e., ML) laser-ranging system. The critical assumption 
made by this program is that the time axis can be divided into a larger number of independent 


*F. Davidson and L. Stephens, “Experimental Performance of Point Process Estimators of Optical Pulse Delay,” to be 
published in IEEE Transactions in Communications, August 1978. 
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time bins, with the average photoelectron emission rate constant over any bin width. As 
the time bins shrink to zero width, the simulation results should approach the analytical 
results. Note that this assumption is consistent with that used in the section on “Mathe- 
matical Representation of a Poisson Process.” The use of a finite number of time bins also 
corresponds to a receiver operating with a finite resolution time and, hence, finite bandwidth. 
Because this will always be the case in practice, this approximation was not believed to be 
restrictive. 

The computer program first generates and stores the optical pulse shape of total energy, 

Q, at each time-bin value. Next, the matched-filter impulse response is generated and stored. 
Then, by considering each time bin independently, a random Poisson process was generated, 
with the mean value of the Poisson process at each time bin equal to the instantaneous 
optical intensity at the same bin. The random-count record was next convolved with the 
matched-filter impulse response to generate the log-likelihood expression of equation 42. 
Finally, the time estimate, r, was chosen as the time bin at which the filter output was 
maximum. This process was repeated many times, while building a histogram of the time 
estimates, r. Finally, the mean and standard deviations of the histogram were calculated. 

The mean value of the histogram, which represented the bias of the simulated ML estimator, 
was always near zero. The standard deviation represented the rms timing error for an ideal 
receiver operating at a given average photoelectron level, Q. By repeating the simulation 
program for different values of Q, the performance of the simulated receiver can be plotted 
as a function of average received photoelectrons per pulse. For programming simplicity, 
only background free cases were studied. 

Figure 9 shows a representative plot of this process for a Gaussian pulse shape and an average 
signal level of 5 photoelectrons per pulse. The time resolution is 40 bins per pulse width and 
has the same scale for all plots. The count record and correlation function shown are those 
for the last (500th) count record, and the histogram is that for the stop times for the random 
Poisson process. The left most histogram bin is that for when no photoelectrons were gen- 
erated in the random count record, and was not used in the calculation of the histogram 
moments. 

Because a stop time was not estimated when no photoevents are recorded, the program com- 
putes the standard deviation of the timing error, given at least one photoelectron occurrence. 
As the width of the time bins shrink to zero, the simulation error results should equal those 
given in equation 65. For large Q, the simulation error results should also equal equation 
61, but discrepancies will exist for small Q values because the approximations used to ob- 
tain equation 61 preclude the case of small Q. 

For program validation, the average number of photocounts after each simulation run was 
compared to the Q value used in the program. In each case, the average number computed 
by the program differed by less than 1 percent from Q. As a second check, the simulated 
performance of a system using the raised cosine pulse defined in equation 62 was compared 
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RECEIVER COUNT RECORD 



HISTOGRAM OE RECEIVER STOP TIMES 



Figure 9c Representative simulation results for Q = 5 and Gaussian pulse shapec 


to the analytical results predicted by equation 65 in the range from 1 to 100 photoelectrons. 
Shown in figure 1 0, these results indicate excellent agreement between analytical and simu- 
lated receiver performance. For all further simulations, a time resolution of 40 bins per 
pulse duration and 500 counts per histogram were also used. 

Figure 1 1 shows simulation results for the ML receiver operating with a Gaussian pulse of 
standard deviation, Up = D/6. For comparison, the simulated performance of a raised cosine 
with a full width of D is also showm Figure 1 1 illustrates the nearly equal timing resolution 
obtainable with these two bell-shaped pulses* 
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Figure 10. Comparison of analytical and simulated ML performance for raised-cosine optical-pulse 
shape, given at least one photoelectron observation and no background counts. 


To investigate the ML receiver performance with rectangular pulse shapes, the raised-cosine 
pulse shape was raised to a fractional power, F. In the limit, as F approaches zero, the pulse 
shape approaches a rectangle. Further simulations were carried out for F values of 0.5, 0.1, 
and OoOl. Figure 12 shows the resulting optical-pulse shapes, and figure 13 shows the per- 
formance of the receiver. These curves, which represent new results in this area, show how 
the slopes of the rms time-jitter curve changes from for bell-shaped pulses, to 

ct^/Q for rectangular pulses, where and are constants. The curve for F = 0.01 is in 
good agreement with 


\/Var{ 

D ^ D 


, V2 


vl 


(1 -6*5-0 6-^5 -Q2e-Q/2) 
2 ( 1 - e-*5 ) 


( 66 ) 


which is the expression calculated by Bar-David for the MMSE timing error when using a 
rectangular optical pulse (Reference 4). 
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Figure 1 1 . Comparison of simulated ML performance for raised-cosine and Gaussian optical-pulse 
shapes, given at least one photoelectron observation and no background counts. 


It is evident from figure 1 3 that, to minimize ranging error, a bell-shaped optical pulse 
shape should be chosen for operating with average signals of less than 15 photoelectrons 
per pulse and that a rectangular pulse should be chosen for operating at higher signal levels. 
The more rapid decrease of ranging error versus received photoelectron level should strongly 
motivate the use of sharp leading- and trailing-edged pulses when possible in pulsed laser 
ranging systems. 

PERFORMANCE OF SUBOPTIMUM RECEIVER STRATEGIES 

Several studies have been made of suboptimal ranging-receiver configurations, including 
systems using a high-resolution leading-edge receiver (References 6 and 7). One stated 
advantage of this and other proposed leading-edge strategies is that they are simplier to 
implement than the ML receiver. Another possible advantage is that electronics are cur- 
rently suited to their implementation or that the strategies can be used when the optical- 
pulse shape is not known in advance, as is the case for satellite-induced signal broadening 
(Reference 12). In this section, the performance characteristics of two possible leading- 
edge receivers are discussed. 
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Figure 12. Pulse shapes used to investigate ML-receiver performance 
changes from bell-shaped to rectangular optical pulses. 


Figure 14 shows the performance of a high time-resolution leading-edge receiver when the 
receiver is operating with a 400-ps rectangular optical pulse. For comparison, figure 14 also 
shows the simulated performance of an ML receiver operating with a nearly identical pulse 
shape (F = 0.01 curve of figure 12). The receiver model for the leading-edge receiver assumed 
a 5-ps time resolution with the trigger threshold set at the 2-photoelectron level. Such a 
receiver can be implemented by doing numerical calculations on the digitized output data 
from a circularly swept image-converter tube. Two such image-converter receiver systems 
are being developed at the Goddard Space Flight Center. 

The performance of this receiver was also analyzed assuming nonideal image-converter gains. 

A gamma distribution was assumed for the pulse-height distribution of the swept-image 
tube, and the receiver performance was calculated for different values of the normalized 
variance of the gamma distribution. Figure 14 shows two performance curves for the 
normalized variance values of the gamma distribution of 1 and 1/3. In the signal region 
shown, the nonideal pulse-height resolution is shown to increase receiver errors only slightly. 
As is readily apparent in the figure, use of such a leading-edge receiver should be justified 
only under the condition of large return-pulse shape uncertainties. When the pulse shape is 
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Figure 13. ML-receiver performance for increasingly rectangular optical-pulse shapes. 


known in advance, such a leading-edge receiver requires approximately 10 times more return- 
pulse energy than an ML receiver to achieve the same ranging performance. 

Figure 1 5 shows the performance of several matched-filter type receivers for a raised-cosine 
optical-pulse shape. The mismatched ML receiver represents a filter matched to the pulse 
shape rather than to the logarithm of the pulse shape. The leading edge receivers have filters 
followed by a discriminator with triggers on the 50-percent risetime point of the filter output 
waveform. The advantage of this type receiver is that such discriminators have been in use for 
many years in both nuclear research and are used in current laser ranging systems (Reference 1). 
The use of such discriminators following a matched-filter degrades receiver error performance 
only sli^tly in the signal region shown. 

The similar performance of the matched-filter and mismatched ML receivers in figure 15 
suggest that this type of receiver structure is relatively insensitive to small optical-pulse 
shape variation. Furthermore, the similar error performance of 50 percent versus peak 
detection also suggests that the exact trigger point as a fraction of the total pulse height 
is not critical if the trigger fraction does not change from one received laser pulse to the 
next. 
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Figure 14. ML and high-resolution leading-edge receiver performance for400-ps rectangular 
optical-pulse shapes. Leading-edge deterministic represents deterministic receiver gain, and K 
is the normalized variance of stochastic receiver gains. 


However, both leading-edge receivers have a signal-energy-dependent time bias, that was not 
present when using peak detection following the matched-filter. This bias consists of a con- 
stant bias that is nominally equal to one-fourth of the filter impulse response full width plus 
a small signal-dependent term. Figure 16 is a plot of this small term. Although the constant 
term of the bias can be removed at any time during operation of the ranging receiver, the 
signal-dependent term could cause significant errors when the desired range resolution is a 
small fraction of the optical-pulse width. By measuring the total number of electrons re- 
ceived per pulse, this smaller signal-dependent bias term can be removed on a pulse-by-pulse 
basis, but only at the expense of additional receiver complexity. Although this receiver 
configuration ig desirable from the viewpoint of range-error performance, these potential 
bias errors might limit the use of constant-fraction discriminators following the receiver 
filter. 
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Figure 15. Performance of several matched-filter type receivers to bell-shaped optical 
pulses, given at least one photoelectron observation^ Mismatched filters represent filter im- 
pulse responses equal to the laser pulse shape rather than to the logarithm of the pulse shape. 
LE represents timing to the 50-percent point on the leading edge of the filter output wave- 
form; the ML receivers time to the peak of the waveform^ 


CONCLUSIONS AND AREAS FOR FUTURE STUDY 

This document has summarized previously proposed methods for optimum direct-detection 
laser-ranging receivers and has compared the performance of these candidate receiver imple- 
mentations by using computer simulations of the detection process. The ML and MMSE 
estimators were found to give the same performance for both bell-shaped and rectangular 
optical-pulse shapes when no background radiation existed. Of several hardware receiver 
implementations presented, the matched-filter peak-detector was found to be preferable 
from the standpoint of optimum performance and receiver simplicity. Although electronics 
are available for using a constant-fraction discriminator following the matched-filter, a 
signal-level-dependent time bias exists with this implementation. To achieve accuracies to 
a small fraction of the optical-pulse width, this bias, which can be compensated for on a 
pulse-by-pulse basis, must be addressed in the receiver design. 
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Figure 16. Signal-dependent range bias when using a matched filter followed by a 50- 
percent risetime discriminator. For symmetrical optical-pulse shapes, the static range bias is 
nominally 1/4D. Note that, for large numbers of photoelectrons per pulse, the range bias can 
be large compared to the random range errors about this value. 


Several areas should be investigated further to fully characterize the ranging performance 
of the receiver configurations addressed in this study. Throughout this document, the 
receiver processing electronics were assumed to have access to the individual photo-occurrence 
times. In turn, these occurrence times were assumed to be unit amplitude delta functions in 
time. In reality, the outputs of the best available optical detectors have both stochastic gains 
and stochastic transit times. These random factors will further degrade the receiver perfor- 
mance from those reported in this document. A more elaborate detection model could take 
these additional receiver uncertainties into account. Further computer simulation work might 
also yield insight into the operation of matched-filter type receivers under nonzero back- 
ground conditions and when using constant-fraction discriminators with sharp-edged optical 
pulses. Detection and falso-alarm probabilities, which have not been addressed here, must 
be studied for implementable receiver strategies. 

Supporting experimental data, especially that with carefully measured average signal levels, 
would be valuable to laser-ranging system designers. Using these data with the analytical 
and simulation studies would not only serve to quantify the range error attributable to each 
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ranging receiver component, but also would lay the required groundwork for much improved 
laser-ranging system performance. 
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